Chapter 5 

Harmonic Oscillator and Coherent 
States 


5.1 Harmonic Oscillator 

In this chapter we will study the features of one of the most important potentials in 
physics, it’s the harmonic oscillator potential which is included now in the Hamiltonian 

V{x) = ^a; 2 . (5.1) 

There are two possible ways to solve the corresponding time independent Schrodinger 
equation, the algebraic method, which will lead us to new important concepts, and the 
analytic method, which is the straightforward solving of a differential equation. 


5.1.1 Algebraic Method 


We start again by using the time independent Schrodinger equation, into which we insert 
the Hamiltonian containing the harmonic oscillator potential (5.1) 


x 2 ) 0 = E-0 . 


. h 2 d 2 mu o 

H * = + ^ " 


We rewrite Eq. (5.2) by defining the new operator x := mux 

2 


(5.2) 


H V = 

2 m 


h d \ , 

—— + (mux) 

i dx 1 


-0 = — \p 2 + x 2 ] -0 = E'lb . 
2m L J 


(5.3) 


We will now try to express this equation as the square of some (yet unknown) operator 

-> (x + ip)(x — ip) = p 2 + x 2 + i(px — xp ), (5.4) 


2 I -2 

p + X 


but since x and p do not commute (remember Theorem 2.3), we only will succeed by taking 
the x — p commutator into account. Eq. (5.4) suggests to factorize our Hamiltonian by 
defining new operators a and cd as: 
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Definition 5.1 


a 


s( mra + P 


annihilation operator 


at := " O 


creation operator 


These operators each create/annihilate a quantum of energy E = hu>, a property 
which gives them their respective names and which we will formalize and prove later on. 
For now we note that position and momentum operators are expressed by a’s and cd’s 

like 

x= \[^( a + a ^ p = ~ i ][^ ( a “ at ) ' ( 5 ‘ 5 ) 


Let’s next calculate the commutator of the creation and annihilation operators, 
quite obvious that they commute with themselves 

[a , a] = [ad , a ' ] = 0 . 

To find the commutator of a with ad we first calculate aaJ as well as a'a 

ad = -- \(mcox) 2 + imu[p, x] + p 2 1 

2 mwn L J 

a'a = -- \(rriLUx) 2 — imw\p , x] + p 2 ] . 

zrmvn L J 


It’s 

(5.6) 


(5.7) 

(5.8) 


Since we know [p, x] = — ih we easily get 

[ a , a ' ] = 1. 


(5.9) 


Considering again the Hamiltonian from Eq. (5.3) we use expressions (5.7) and (5.8) 
to rewrite it as 

H = —*— (p 2 + (rmvx) 2 ) = — (a)a + acd) . (5.10) 

2 m 2 v ' 

LIsing the commutator (5.9) we can further simplify the Hamiltonian 

[a, ad] = a a) — a) a = 1 a a) = a) a + 1 , (5-11) 


H = hj{a)a + |) . (5.12) 

For the Schrodinger energy eigenvalue equation we then get 

huj{a)a + |) i/j — Eip, 


(5.13) 
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which we rewrite as an eigenvalue equation for the operator a'a 


a f o. L = I --- 1 ip , 


E 1 


(5.14) 


, fruj 2 

having the following interpretation: 

Definition 5.2 N := a) a occupation (or particle) number operator 


and which satisfies the commutation relations 

[N , aJ ] = a) [IV, a] = —a 


(5.15) 


Next we are looking for the eigenvalues v and eigenfunctions if) v of the occupation 
number operator N, i.e. we are seeking the solutions of equation 


N ip v = uip v . 


(5.16) 


To proceed we form the scalar product with ip u on both sides of Eq. (5.16), use the positive 
definiteness of the scalar product (Eq. (2.32)) and the definition of the adjoint operator 
(Definition 2.5) 


v ( tj) v | — (i>„\ N | a f a ) = (a \ aij) v ) > 0 . 


(5.17) 


The above inequality vanishes only if the corresponding vector a ip is equal to zero. Since 
we have for the eigenvalues v > 0, the lowest possible eigenstate corresponds to the 
eigenvalue v — 0 

v — 0 —» o = 0. (5.18) 

Inserting the definition of the annihilation operator (Definition 5.1) into condition 
(5.18), i.e. that the ground state is annihilated by the operator a, yields a differential 
equation for the ground state of the harmonic oscillator 


o 


1 

\/ 2 mcJh 


. h d . , 

mux + i -) 

i ax 


=> 


mu 


h 


-x + 



^0 


0 

0 . 


(5.19) 


We can solve this equation by separation of variables 


‘ dil) o 

% 


= — dx 


mu 


h 


x 


, , mu 2 ^ ^r 

hi '0o = - x + In J\f , 
2 n 


(5.20) 


where we have written the integration constant as In A/", which we will fix by the normal¬ 
ization condition 
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l "(^) = ~W x2 * MX) • ( 5 - 21 ) 

We see that the ground state of the harmonic oscillator is a Gaussian distribution. The 

OO 

normalization f dx \'i/jo(x)\ 2 = 1 together with formula (2.119) for Gaussian functions 

— OO 

determines the normalization constant 


u 2 = 



u 



1 

4 


(5.22) 


We will now give a description of the whole set of eigenfunctions ij) v of the operator 
N based on the action of the creation operator by using the following lemma: 


Lemma 5.1 If if v is an eigenfunction of N with eigenvalue v, then 
of il) u also is an eigenfunction of N with eigenvalue (y + 1). 


Proof: 


Na^„ Eq '= 15) (a'N + af)^ v = of {N + l)if„ 

= a' {u + 1) — (y + 1) a* if v . q.e.d. 


(5.23) 


The state a'rf v is not yet normalized, which means it is only proportional to if v + \ , 
to find the proportionality constant we use once more the normalization condition 

( | of'f’v ) = {ipu\ acf if v ) = {if u \ (N + 1) | if v ) = (u + 1) ( if v \ if v ) (5.24) 


at a + 1 


a 


= Vi; + l Vwi • 


(5.25) 


This means that we can get any excited state of the harmonic oscillator by succes¬ 
sively applying creation operators. 


In total analogy to Lemma 5.1 we can formulate the following lemma: 


Lemma 5.2 If if v is an eigenfunction of N with eigenvalue u, then 
aif v also is an eigenfunction of N with eigenvalue [v — 1). 
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Proof: 


N av u Eq ' - l5) (aN — a)ip u = a(N — 1) v„ 

= a (u — 1) — (y — 1) a'lpu . q.e.d. 


We again get the proportionality constant from the normalization 


(5.26) 


( aif u | aif> v ) = ( f/y | oTo, if v ) = (ip„\ N \ij; u ) = v^ v \^ v ) 

N v\i>0 1 


(5.27) 


a'ip v = \/v ipu-i . 


(5.28) 


Summary: Eigenvalue equation for the harmonic oscillator 

The time independent Schrodinger equation, the energy eigenvalue equation 

H^ n = E n i> n , (5.29) 

with the Hamiltonian H = hu (.N + |) and the occupation number operator N = cda 
provides the energy eigenvalues 


E n = fmj {n + \) . (5.30) 

The corresponding state vectors, the energy eigenfunctions, are given by 

Mx) = 7^ {atrMx) = 7^ S) 4 (at) ° eM ~'fK x2) ’ (5 ' 31) 

with the ground state 

i 

M x ) = (^r) ) • (5-32) 

The explicit form of the excited state wave functions will be calculated later on but we 
can for now reveal that they are proportional to a product of the ground state and a 
family of functions, the so-called Hermite polynomials H n . 

The wave functions thus form a ladder of alternating even and odd energy states, see 
Fig. 5.1, which are each separated by a quantum of energy hu, i.e. equally spaced. The 
creation and annihilation operators then ’’climb” or ’’descend” this energy ladder step by 
step, which is why they are also called ladder operators. 



100 


CHAPTER 5. HARMONIC OSCILLATOR AND COHERENT STATES 



Figure 5.1: Harmonic oscillator: The possible energy states of the harmonic oscillator 
potential V form a ladder of even and odd wave functions with energy differences of 

huo. The ground state is a Gaussian distribution with width x 0 = \J picture from 
http: //en. wikipedia.org/wiki/Quantum_mechanical_harmonic_oscillator 


5.1.2 Zero Point Energy 

We already learned that the lowest possible energy level of the harmonic oscillator is not, 
as classically expected, zero but E 0 = \frw. 

It can be understood in the following way. The ground state is an eigenfunction of 
the Hamiltonian, containing both kinetic and potential energy contributions, therefore 
the particle has some kinetic energy in the vicinity of x — 0, where the potential en¬ 
ergy V{x — > 0) —> 0. But this implies according to Heisenberg’s uncertainty relation 
(Eq. (2.80)) that the momentum uncertainty increases pushing the total energy up again 
until it stabilizes 

A p oc /* for Ax —> 0 ^ / 0. (5.33) 

Ax 

We will now illustrate the harmonic oscillator states, especially the ground state and 
the zero point energy in the light of the uncertainty principle. We start by calculating 
the position and momentum uncertainties using Definition 2.10 

Eq (5 5) 

( x ) = ( V’n | X | V’n ) 'oc ' ( ljj n | a + a t I ) OC ( V> n I V’n-l ) + ( i>n | V’n+l ) = 0 

"-v-' "•-v-' 

0 0 

(5.34) 

h 

(x 2 ) = (-0n| X 2 |-0n) = -- ('Ipn \ 0? + dd) + Cl)CL + (cJ) 2 | ) . (5.35) 

Imuj 

Because the the eigenstates for different n are orthogonal, the expectation values 
of a 2 and (aQ 2 vanish identically and we proceed by using Eq. (5.11), where a)a = N. 


















5.1. HARMONIC OSCILLATOR 


101 


We get 


(x 2 ) = ( ip n | 27V + 1 | ip n ) = —(n + l) = x 2 0 (■n + \) , 

x / 2mu mu z 


(5.36) 


where we have introduced a characteristic length of the harmonic oscillator = y ~ 
and n — 0,1, 2, • • • is a natural number. The position uncertainty is then given by 


(Aa;) 2 = ( x 2 ) - ( x ) 2 = x 2 0 (n + \) , 
which vanishes in the classical limit, i.e. h —> 0 or m —> oo. 


(5.37) 


Let us continue with the momentum uncertainty 


Eq (5 5) 

(p) = (V’n | P | V’n) OC' (V’nl O - O 1 * |^n) OC (V’nl V’n-l) - (V’nl f/'n+l) = 0 


(5.38) 


(/) = (V’n | P 2 |^n) = 


ti 2 , . 

— W ( n + 2) ’ 


muh 


(ip n | a" — am — da + (a') | ip n ) (5.39) 

-(2JV+1) 


where we have again used the characteristic length xq, the ladder operator commuta¬ 
tion relations (Eq. (5.11)) and the orthogonality of the eigenstates. For the momentum 
uncertainty we then have 


(A p) 2 = ( p 2 ) - ( p ) 2 = -2 (n + |) • (5.40) 

x o 

Constructing the uncertainty relation from Eq. (5.37) and Eq. (5.40) we finally get 


Ax Ap — h(n + |) > | generally (5.41) 

= | for the ground state n — 0 . (5.42) 

The ground state, Eq. (5.32), can be rewritten in terms of the characteristic length Xq 

1 x 2 

M x ) = 7-= exp(--^) , (5.43) 

VA) V 71- za; o 

which is a gaussian wave packet with standard deviation xq characterizing the position 
uncertainty in the ground state, see Fig. 5.1. 

Since we have made clear that the zero point energy is in accordance with the un¬ 
certainty relation we can even strengthen this correlation in so far as the non-vanishing 
zero point energy is a direct consequence of the uncertainty principle. Starting from the 
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uncertainty relation (Eq. (2.80)) and keeping in mind that here the expectation values of 
the position- and momentum operator vanish identically, see Eq. (5.34) and Eq. (5.38), 
we can rewrite the uncertainty relation as 


(x 2 ) (p 2 ) > \ . 


(5.44) 


We now calculate the mean energy of the harmonic oscillator, which is the expectation 
value of the Hamiltonian from Eq. (5.2) 


E = {H) — 


( P 2 ) 

2m 


moo~ 


2 v Ec »- 0 5 - 44 ) ( p 2 ) i mu 2 K 2 1 


x ) > 


2m 


+ —— -r 


4 {p 2 ) 


(5.45) 


To End the minimal energy we calculate the variation of Eq. (5.45) with respect to ( p 2 ), 
which we then set equal to zero to hnd the extremal values 


1 

2 m 


moo 


2 h 2 l 


= o =* 


p 2 ) . = 

1 / mm 


mooh 


2 4 (p2) 2 ' XJr 7 mm 2 

Inserting the result of Eq. (5.46) into the mean energy (Eq. (5.45)) we get 

2 hjj hoo hoo 


(5.46) 


rn 1 mooh moo 

E > -+ - 

~ 2m 2 8 


2 ft 2 


mooh 4 4 2 


thus 


E> 

~ 2 


(5.47) 

(5.48) 


Theorem 5.1 (Zero point energy) 

The zero point energy is the smallest possible energy a physical system can 
possess, that is consistent with the uncertainty relation. It is the energy of its 
ground state. 


5.1.3 Comparison with the Classical Oscillator 

In this section we want to compare the quantum oscillator with predictions we would get 
from classical physics. There the classical motion is governed by the position function 


x(t) = qo sin(caf), 

and the total energy of the system is constant and given by 


(5.49) 


E = 


moo 


% 


2 


(5.50) 
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where in both cases, the angular frequency cu is fixed and g 0 is a given starting value. 
To make a comparison with quantum mechanics we can construct a classical probability 
density l / i' 7 c iass(^') in the following way 

W c i ass (x)dx = (5.51) 

where dt is the time period where the classical oscillator can be found in the volume 
element dx, and T = — is the full oscillation period. We then calculate the volume 
element dx by differentiating equation (5.49) 

dx = qoiu cos (cut) dt = qoco yjl — sin 2 (cut) dt Eq ' if- 49 ) q QU1 ^J\ — (|-) 2 dt . (5.52) 

Finally inserting Eq. (5.52) into Eq. (5.51) we get the classical probability density of 
the oscillator as a function of x, see Fig. 5.2 (dashed curve) 

^dass(x) = (2vrg 0 ^l - (f) 2 )" 1 • (5.53) 



Figure 5.2: Oscillator probabilities: Comparison of the quantum probability (solid curve), 
\ \ 2 in case of n = 100, with a classical oscillator probability (dashed curve), Eq. (5.53). 
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The starting value g 0 is fixed by the comparison of the energy of the quantum oscillator 
with the corresponding classical energy 

„ 2n +1 mu 2 2 

En = - - - ftw < -> Wlass = ~Y~ % ■ ( 5 - 54 ) 

It allows us to relate the starting value qo of the classical oscillator to the characteristic 
length Xo of the quantum mechanical oscillator 

(2 n + 1 )h , - 

--— = y/2n + 1 x 0 . 5.55 

mu 

In the comparison quantum versus classical oscillator we find the following. For the 
low-lying states, of course, the two probabilities for finding the particle differ considerably 
but for the high-lying states the smoothed quantum probability resembles very much the 
classical one, as can be seen in Fig. 5.2 for the case n = 100. 



5.1.4 Matrix Representation of the Harmonic Oscillator 

In this section we now want to briefly sketch how the harmonic oscillator problem can 
be written in a matrix formulation, which introduces us to the concept of the Fock space 
and the occupation number representation. It is the appropriate formalism for relativistic 
quantum mechanics, i.e. Quantum Field Theory (QFT). 

Starting from the eigenfunctions of the harmonic oscillator {^ n (a;)}, which form a 
complete orthonormal system of the corresponding Hilbert space, we change our notation 
and label the states only by their index number n, getting the following vectors 1 


4>o 


fl\ 


fo\ 


: 

0 

i i 

1 


0 

0 

T 

O 

0 

= 1 1 ) , -»■ 

1 

• 

1 1 


/ 

0 

\7 

\7 

n— th row 

i : i 




n) . (5.56) 


The vectors { | n) } now also form a complete orthonormal system, but one of the so 
called Fock space or occupation number space 2 . We will not give a mathematical definition 
of this space but will just note here, that we can apply our rules and calculations as before 
on the Hilbert space. The ladder operators however have a special role in this space as 
they allow us to construct any vector from the ground state 3 , see Eq. (5.58). 


technically the vectors of Eq. (5.56) are infinite-dimensional, which means that the column-matrix 
notation is problematic, the ket-notation however is fine and we accept the column-vectors as visualization. 

2 In many body quantum mechanics and QFT the occupation number n is no mere label for the physical 
states, but describes the particle number of the considered state. This means that in QFT particles are 
described as excitations of a field, composed of a multitude of harmonic oscillators. The creation- and 
annihilation operators thus create and annihilate particles. 

Which in QFT is the vacuum state, often also denoted by | O). 
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Let us state the most important relations here: 

( n\m) = S n>m 


n) = 


(o'rio) 


a' I n 


) = y/n + 1 | n + 1) 
a | n) = yfn \ n — 1) 

N | n ) = a'a | n ) = n \ n) 

H \ n) = Huj (N + |) | n ) = ftw (n + |) | n ) . 


(5.57) 

(5.58) 

(5.59) 

(5.60) 

(5.61) 

(5.62) 


With this knowledge we can write down matrix elements as transition amplitudes. We 
could, for example, consider Eq. (5.57) as the matrix element from the n-th row and the 
m-th column 


(5.63) 



(l 0 • 


n \ m) = 5 n)in <—> 

0 1 



V 

•7 


In total analogy we can get the matrix representation of the creation operator 


/ 0 0 0 ••• \ 
^10 0 
0 V2 0 
: Vs 

\ / 


(k | a 1 ' | n) — \fn + 1 (k \ n + 1) *—> a 1 ' = 

<5k,n+l 


(5.64) 


as well as of the annihilation operator, occupation number operator and Hamiltonian 

•• \ 



/o v 7 ! o ••• 
0 0 V2 
0 0 0 


V 


(5.65) 


( k | N | n) — n (k \ n) *—> N = 


/0 0 0 
0 1 0 
0 0 2 


\ 


(5.66) 
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{k | H | n) = hw(n + |) (k\n) 


I 0 0 

0 I 0 
o o I 


. (5.67) 


5.1.5 Analytic Method 

In this last section on the harmonic oscillator we want to explore the straightforward 
method of solving the Schrodinger equation as a differential equation 




(5.68) 


Multiplying both sides by ^ and introducing the dimensionless variables £ = ^ and 
K = where we again have = \/—, we obtain the differential equation 


= (e - k) m ■ 


(5.69) 


To get an idea how the solution could look like, we first try to solve this equation for 
large £, which means we can neglect the K -terms 


^£2 = £ 2 ^(0 • 

For this differential equation we try the ansatz i/j = exp(—y), providing 

= -? exp(-k) , A vk) = (y - i) vk) • 


(5.70) 


(5.71) 


Since we are in the limit of large £, we can drop 1 in the last term on the right hand side 
of Eq. (5.71) and we have made the correct ansatz (in this approximation). To End the 
general solution we can use a similar ansatz, differing only insofar as we multiply it with 
a function of £, a method known as variation of the constants 


m = exp(-V) . 


Differentiating our ansatz (5.72) we obtain 


= 

d 2 

- 


T -«M0)exp(-|) 


(5.72) 


(5.73) 


d£2 


dh(£) 


+ (e 2 - l)h(0 exp(-V). (5.74) 
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On the other hand we know the second derivative of if from the Schrodinger equation 
(5.69), leading to a differential equation for the function h(£) 


d 2 h(f) 

d£ 2 


- + (K - l)h(() 


0 . 


(5.75) 


which can be solved by a power series ansatz. We are not computing this here, but direct 
the interested reader to a lecture and/or book about differential equations 4 . Equation 
(5.75) is the differential equation for the so called Hermite polynomials H n , thus we 
replace our function by 


MO 


H n (0 , n = 0 , 1 , 2 ,--- 


(5.76) 


where (.K — 1) = 2n or K = = 2 n + 1 respectively. With this substitution we directly 

find the energy levels of the harmonic oscillator, as calculated earlier in Eq. (5.30) 

E n — fuv (n + |) . (5-77) 

With the help of the Hermite polynomials we can now also give the exact expression 
for the eigenfunctions of the harmonic oscillator 


l 


i’n = 




rnu \ 4 X x~ 

_j ^(_ )exp( __) 


(5.78) 


[nl) wift/ xo 

The Hermite Polynomials can be easiest calculated with the Formula of Rodrigues 

H a (x) = (-ire- 2 . (5.79) 

To give some examples we list some of the first Hermite Polynomials 

Hq{x) = 1, Hi(x) = 2x , H 2 (x ) = 4x 2 — 2 

H 3 (x) = 8x 3 - 12.r| H 4 (x) = 16a; 4 - 48a; 2 + 12 . (5.80) 

A last interesting property is the orthogonality relation of the Hermite Polynomials 


dx H m (x) H n (x) e x = y/n 2 n (n!) 8 


(5.81) 


5.2 Coherent States 

Coherent states play an important role in quantum optics, especially in laser physics and 
much work was performed in this field by Roy J. Glauber who was awarded the 2005 
Nobel prize for his contribution to the quantum theory of optical coherence. We will try 
here to give a good overview of coherent states of laser beams. The state describing a 
laser beam can be briefly characterized as having 


4 See for example [14]. 
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1 . an indefinite number of photons, 

2 . but a precisely defined phase, 

in contrast to a state with fixed particle number, where the phase is completely random. 

There also exists an uncertainty relation describing this contrast, which we will plainly 
state here but won’t prove. It can be formulated for the uncertainties of amplitude and 
phase of the state, where the inequality reaches a minimum for coherent states, or, as we 
will do here, for the occupation number N and the phase 5 

AN A(sin<3>) > - cos<f> , (5.82) 

which, for small <f>, reduces to 

1 

AIV A<f> > - . (5.83) 

5.2.1 Definition and Properties of Coherent States 

Since laser light has a well-defined amplitude - in contrast to thermal light which is a 
statistical mixture of photons - we will define coherent states as follows: 

Definition 5.3 A coherent state \ a ) , also called Glauber state, 

is defined as eigenstate of the amplitude operator, 
the annihilation operator a, with eigenvalues a € C 

a | a) = a \ a) . 


Since a is a non-hermitian operator the phase a = |a| e llf € C is a complex number 
and corresponds to the complex wave amplitude in classical optics. Thus coherent states 
are wave-like states of the electromagnetic oscillator. 


Properties of coherent states: 

Let us study the properties of coherent states. 

Note: The vacuum | 0 ) is a coherent state with a = 0 . 

Mean energy: 

( H ) = ( a | H | a) = hw (a \ a^a + \ \ a ) = hw (\a\ 2 + |). (5.84) 

The first term on the right hand side represents the classical wave intensity and the second 
the vacuum energy. 


5 See for example [15]. 
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Next we introduce the phase shifting operator 

U{9) = e~ i9N , (5.85) 

where N represents the occupation number operator (see Definition 5.2). Then we have 
U\6)aU{6) = ae~ ie or U{6) a U\6) = a e ie , (5.86) 

i.e. U ( 6 ) gives the amplitude operator a phase shift 9. 

Proof: The l.h.s. ±U*{9)aU{9) = iW(6) [N,a]U(9) = -i U\9) a U{6) and the 

r.h.s. -jnae~ t6 = —i a e~ ld obey the same differential equation. 


The phase shifting operator also shifts the phase of the coherent state 

U{9) | a) = | ae~ id ) . (5.87) 

Proof: Multiplying Definition 5.3 with the operator U from the left and inserting the 

Unitarian relation U'U — 1 on the l.h.s we get 

IJaU^ U \a) = aU \ a) 

ae‘ # 

=>• aU \ a) = ae~' ld U |a) . (5.88) 

Denoting U \ a) = | a ) we find from Definition 5.3 that a | a ) = a | a ) and thus 
a = ae , q.e.d. 


Now we want to study the coherent states more accurately. In order to do so we 
introduce the so called displacement operator which generates the coherent states, similar 
to how the creation operator a7, see Eq. (5.58), creates the occupation number states | n). 

Definition 5.4 The displacement operator D(a) is defined by 

D(a ) = e aa) ~ a * a 

where a = |a| e lif G C is a complex number and af 
and a are the creation - and annihilation - operators. 


The displacement operator is an unitary operator, i.e. D^D = 1 and we can also 
rewrite it using Eq. (2.72). To use this equation we have to ensure that the commutators 
[[ A , B] , A] and [[A, B] , B] vanish, where A = a of and B = a* a . We therefore 
start by calculating the commutator of A and B 

[A, B] = [ a a' y a* a ] = a a* [ a' , a ] = — \a 


-l 


(5.89) 
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Since the result is a real number, it commutes with A and B and we are allowed to use 
Eq. (2.72) where we can immediately insert the result of Eq. (5.89) to get the displacement 


operator in the following form: 

D(a) = e"^ |a|2 e aat e a * a . (5.90) 

Properties of the displacement operator: 

I) D\a) = D~ 1 (a) = D(—a) unitarity (5.91) 

H) D\a) a D(a) = a + a (5.92) 

HI) D ] {a)a) D{a) = + a* (5.93) 

IV) D(a + (3) = D{a) D((5) e~ ilm(a ^ (5.94) 


We will prove two of these statements and leave the other two as an exercise. 


Proof: of II) 

D\a)aD{a) = e a * a ~ aa] ae aat ~ a * a = a + [a*a — ao^, a] = 

= a + a* [a , a] — a [a' , a] = a + a q.e.d. (5.95) 

o -1 

We have used the Baker-Campbell-Hausdorff formula (Eq. (2.71)), where A = a*a — a a' 
and B = a. The higher order commutators vanish, since the commutator of A and B is 
a complex number, that commutes with the other operators. 


Proof: of IV) 

D(a + (5) = 


aa^ — a*a-\-P —/3*a aa^ — a*a (3 —/3*a — 7 ; \ — a*a,/3 —/3*a] 

c- — C- C- c- ^ L -* — 

D{a)D(f5)e-\ {a ^- a ^ ] = D(o)D(f5)e- ilra{a ^ ) q.e.d. (5.96) 


Here we have used Eq. (2.72), with A — a cd — a*a and B = (3 cd — /3*a which is justified 
because the condition that the higher commutators vanish is satisfied since [A , B] e C. 
The result of this last commutator can also be seen easily by noting that the creation- 
and annihilation-operators commute with themselves but not with each other, giving ±1. 



5.2. COHERENT STATES 


111 


With this knowledge about the displacement operator we can now create coherent 
states. 


Theorem 5.2 The coherent state \ a) is generated from the vacuum | 0 ) by 
the displacement operator D(a ) 

| a ) = D(a ) | 0} . 


The ” vacuum” j 0) is the ground state with occupation number n = 0, which is defined 
by a | 0) = 0, see also Eq. (5.18). 


Proof: Applying a negative displacement to | a) we find from properties (5.91) and 

(5.92) that 

a D(—a) | a) = D(—a ) D\—a) a D(—a) \ a ) 

' -V-' 

a — a 

=>- aD(—a ) | a) = D(—a)(a — a) |a) = 0. (5.97) 

The r.h.s. vanishes because of Definition 5.3. This implies that D(—a) \ a) is the vacuum 
state | 0) 

D(—a) | a) = 10) l c 0 = D(a) |0) q.e.d. (5.98) 


On the other hand, we also could use Theorem 5.2 as a definition of the coherent states 
then our Definition 5.3 would follow as a theorem. 


Proof: 


0 ) 

= a D(a ) | 0 ) = a | a ) 

a \a) = a \ a) q.e.d. (5.99) 

The meaning of our Definition 5.3 becomes apparent when describing a coherent pho¬ 
ton beam, e.g. a laser beam. The state of the laser remains unchanged if one photon is 
annihilated, i.e. if it is detected. 



D(a) D\a ) a \ a) = D(a) D\a) a D(a ) | 0 ) = D(a) ( a j 0) + a 
1 




112 


CHAPTER 5. HARMONIC OSCILLATOR AND COHERENT STATES 


Expansion of coherent state in Fock space: The phase | a) describes the wave 

aspect of the coherent state. Next we want to study the particle aspect, the states in Fock 
space. A coherent state contains an indefinite number of photons, which will be evident 
from the expansion of the coherent state | a) into the CONS of the occupation number 
states (|n)}. Let us start by inserting the completeness relation of the occupation 
number states 

| a ) = | n) ( n \ a) . (5.100) 

n 

Then we calculate the transition amplitude (n\a) by using Definition 5.3, where we 
multiply the whole eigenvalue equation with (n | from the left 


(n | a | a ) = (n \ a \ a) . 

Using the adjoint of relation (5.59) we rewrite the left side of Eq. (5.101) 


a' I n 


) = \Jn + 1 | n + 1) 


(n | a — \Jn + 1 (n + 1 


\Jn + 1 (n + 1 | a) = a (n \ a) , 

and replace the occupation number n with n — 1 to obtain 

, | , a . . 

(n | a) = —j= {n — 1 | a ) . 

\ n 


By iterating the last step, i.e. again replacing n with n — 1 in Eq. (5.104) 
the result on the right hand side, we get 


( n | a) = 


a 


V n ( n - !) 


( n — 2 | a) = 


a'" 


' 0 | a) 


which inserted into the expansion of Eq. (5.100) results in 


«) = (° k) 


a 


n =0 


n) . 


(5.101) 

(5.102) 

(5.103) 

(5.104) 
and reinserting 

, (5.105) 

(5.106) 


The remaining transition amplitude (0 | a) can be calculated in two different ways, 
first with the normalization condition and second with the displacement operator. For the 
sake of variety we will use the second method here and leave the other one as an exercise 

(0 | a) = (0| D(a) |0) Eq 'i 5 - 90) ( 0 | e Qat |0) . (5.107) 

We expand the exponentials of operators in the scalar product into their Taylor series 

(0| e aaf e a * a |0) = (0| (1 + aal + • • • )(1 + a" 1 a + • • •) | 0) = 1 , (5.108) 

and we let the left bracket act to the left on the bra and right bracket to the right on the 
ket. In this way we can both times use the vacuum state property a \ 0) = 0 to eliminate 
all terms of the expansion except the first. Therefore the transition amplitude is given by 

(0 1a) = e -iM\ 


(5.109) 
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and we can write down the coherent state in terms of an exact expansion in Fock space 


a) 


e 



OO 


E 

n =0 



n) 


e -h a \ 2 




(5.110) 


Probability distribution of coherent states: We can subsequently analyze the 

probability distribution of the photons in a coherent state, i.e. the probability of detecting 
n photons in a coherent state j a ), which is given by 

o |cd 2n e~\ a \ 

Pin) = \(n\a)\ = ——- . (5.111) 

n! 

By noting, that the mean photon number is determined by the expectation value of 
the particle number operator 

n — (a | N | a ) = (a | ata \ a ) De = 5 ' 3 |a| 2 , (5.112) 

we can rewrite the probability distribution to get 

H p Tl 

P{n) = , (5.113) 

nl 


which is a Poissonian distribution. 

Remarks: First, Eq. (5.112) represents the connection between the mean photon 

number - the particle view - and the complex amplitude squared, the intensity of the 
wave - the wave view. 

Second, classical particles obey the same statistical law, the Poisson formula (5.111), 
when they are taken at random from a pool with |cr| 2 on average. Counting therefore 
the photons in a coherent state, they behave like randomly distributed classical particles, 
which might not be too surprising since coherent states are wave-like. Nevertheless, it’s 
amusing that the photons in a coherent state behave like the raisins in a Gugelhupf 6 , which 
have been distributed from the cook randomly with |a| 2 on average per unit volume. Then 
the probability of finding n raisins per unit volume in the Gugelhupf follows precisely law 
(5.111). 

Scalar product of two coherent states: From Theorem 4.3 we can infer some im¬ 

portant properties of the particle number states, which are the eigenstates of the particle 
number operator N = cbo. Since this operator is hermitian, the eigenvalues n are real and 
the eigenstates | n) are orthogonal. But the coherent states are eigenstates of the annihi¬ 
lation operator a, which is surely not hermitian and as we already know, the eigenvalues 
a are complex numbers. Therefore we cannot automatically assume the coherent states 


6 A traditional Austrian birthday cake. 
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to be orthogonal, but have to calculate their scalar product, using again Theorem 5.2 and 
Eq. (5.90) 

(P | a) = (0| D\(3)D(a) |0) = (0| e~ 0at e 0 * a e aat e - a * a |0) e ~ |<l«l 2 + l^l 2 ) . ( 5 . 114 ) 

We now use the same trick as in Eq. (5.108), i.e., expand the operators in their 
Taylor series, to see that the two ’’outer” operators e _da = (1 — /3cd + • • •) and e~ a * a = 

(1 — a*a H-), acting to left and right respectively, annihilate the vacuum state, save for 

the term ”1” and we can thus ignore them. We then expand the remaining operators as 
well and consider their action on the vacuum 

(0 | (1 + /3*a+ — (/3*a ) 2 + •••)(! d - 01 a) + —(a g ^) 2 + •••) |0) (5.115) 

= (••• + (2 I ^ V 2 \(( 3*) 2 + < 1 | ( 3 * + (0 I) (I 0) + a | 1 > + ^v^a 2 |2) + •••) . 

Using the orthogonality of the particle number states (n \m) = <5 nm we get 
(p | a) = e -| (|a|2 + l / 3 l 2 ) (1 + ap* + ^(a/T ) 2 + •••) = e - i (|a|2 + l/3|2) + a/3 * . (5.116) 
We further simplify the exponent, by noting that 

| a — P \ 2 = (a — P)(a* — P*) = \a\ 2 + \P\ 2 — aP* — a*P , (5.117) 

and we can finally write down the transition probability in a compact way 

\(p\a)\ 2 = e~\ a - 0 \ 2 . (5.118) 

That means, the coherent states indeed are not orthogonal and their transition prob¬ 
ability only vanishes in the limit of large differences | a — P \ 1 . 


Completeness of coherent states: Although the coherent states are not orthog¬ 

onal, it is possible to expand coherent states in terms of a complete set of states. The 
completeness relation for the coherent states reads 


1 

7T 



a)(a | = 1 . 


(5.119) 


In fact, the coherent states are ’’overcomplete”, which means that, as a consequence 
of their nonorthogonality, any coherent state can be expanded in terms of all the other 
coherent states. So the coherent states are not linearly independent 


Jd 2 a \ a) (a \ P) 



a)e-i<H 2 + l/?l 2 ) + ^‘ _ 


\P) 


71 


(5.120) 
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Proof of completeness: 


— d 2 a | a) (a 


Eq. (5.110) 1 


7T 


y 1 

7 T ^ y/n\ ml 


n 


) (m | d 2 a e a n (a*) m . (5.121) 


n,m 


The integral on the right hand side of Eq. (5.121) can be solved using polar coordinates 


a = ae 1 *’ = re t<p 


d 2 a e a n (a*) m = / r dr e r ~ r 


2 n oo 

n-\-m f 


dip e* 


o 


= 2n J r dr e r r 
o 


2?r Snm 


Substituting r 2 = t, 2 rdr = dt to rewrite the integral 


OO OO 

—r 2 ,„2n f 4 n 


2 J rdre r r 2n — J e 1 t n dt = T(n + 1) = n\ , 
o o 


(5.122) 


(5.123) 


we recover the exact definition of the Gamma function, see [13, 14]. With this we can 
finally write down the completeness relation 

\ f \ \ _ 

— / d 2 a I a) (a I = — — \n) (n\ nn\ = V I n) (n I = 1 q.e.d. (5.124) 

7 T J 7T Til 


5.2.2 Coordinate Representation of Coherent States 

Recalling the coherent state expansion from Eq. (5.110) 


v -k\a \ 2 a I \ -^|o | 2 

«) = e 2 ' 1 y — \ n ) = e 2 i I ^ 
vn! 


a a 


ty 


n! 


0 ) , 


(5.125) 


n=0 v n =0 

and the coordinate representation of the harmonic oscillator states (Eq. (5.78)) 


'ijjn(x) = (x \ n) = ip n (x) = 


e 


y / 2 n (n\)y/TTX 0 


H n (0 e-T 


l _ i! 

^o(x) = (x|0) = 'ip 0 (x) = i=^= e 2 , 


(5.126) 

(5.127) 


where £ = ^ and Xo = y 777’ we can eas ily gi ye th e coordinate representation of the 
coherent states 

^ OO / "J" \ 7T, 

(x |o) = = e- 2 '“P V 4=V>nW = e-2'“ p W . (5.128) 

ti ^ n ' 
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Time evolution: Let’s first take a look at the time evolution of the harmonic oscillator 

states and at the end we recall Theorem 4.1 and Definition 4.1. Using the harmonic 
oscillator energy (Eq. (5.77)) we easily get 


i>n(t,x) = ipn(x)e n Ent = ^n{x)e~ 


iujt 

e~ 2 


(5.129) 


We are then in the position to write down the time evolution of the coherent state 


J- I lO L w -- 

<i> a (t,x) = e _ 2l Q l e~ 2 


iu>t (ae~ iuJt ) n 


n —0 


\fn\ 


1pn(x) ■ 


(5.130) 


With a little notational trick, i.e. making the label a time dependent a(t) = ae lwt 
we can bring Eq. (5.130) into a more familiar form 


iutt 

0a(t,x) = 0 a ( t )(z) e_2 > (5.131) 

which we can identify as a solution of the time dependent Schrodinger equation. 


Expectation value of x for coherent states: Subsequently, we want to compare 

the motion of coherent states to that of the quantum mechanical (and classical) harmonic 
oscillator, which we will do by studying the expectation value of the position operator. 
For the quantum mechanical case we have already calculated the mean position, which 
does not oscillate (see Eq. (5.34)) 


( x )o,an,tor = 0- (5-132) 

Recalling the expression of x in terms of creation and annihilation operators, Eq. (5.5), 
and the eigenvalue equation of the annihilation operator, Definition 5.3, we can easily 
compute the expectation value of the position for the coherent states 

) coherent = { ( t ) a( t) | X | 0 a ( t ) ) = — y= ( \ a + a ^ | 4>a( t) ) = (5. 

= —^= ( a(t ) + a*(t)) = V2xoRe(a(t)) = V2x 0 |a| cos(cnt — (p ) , 

V 2 

where we used a(t) = ae - " 1 = |a| To summarize the calculation, we conclude 

that the coherent state, unlike the quantum mechanical harmonic oscillator, does oscillate , 
similar to its classical analogue 



1 ^ = y/2 Xq \a\ cos(o;t 



coherent 


(5.134) 
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Coordinate representation in terms of displacement operator: In this section 
we want to express and study the coordinate representation of coherent states via the 
displacement operator D 

4> a (x) = (x | a) = (x| D(a) | 0) = (x \ e aat ~ a * a j 0) . (5.135) 

To simplify further calculations we now rewrite the creation and annihilation operators 
from Definition 5.1 in terms of the dimensionless variable £ = —, 

° = ^ + d£ } ’ at = ^ ~ d£ } • (5 - 136) 

This rewritten operators (5.136) we insert into the exponent of the displacement operator 
from Eq. (5.135) 


t * a fy d \ °* d 

““ " “ “ = - 7! ({ + = 

— -4= ( a - i'M - ~1=4 :( q + o' j — 

22 lm(a) 2Re(o;) 

= V^iIm(o;)£ — v / 2Re(o;) — , 


(5.137) 


and get for the wavefunction (5.135) the following expression 


( (x) = (x | a ) = e 


\/2 i Im(a) £ — \/2 Re(a) 




(5.138) 


where 'ipo is the harmonic oscillator ground state (Eq. (5.127)). 
The time dependent version follows outright 


(j) a (t,x') 0a(t) 


icot 

e~ 2 = 


iujt 


iujt 


(x\a(t))e~r =e-T<z| e Q (*) aW (*) a |0) 




e -^ e ^*Im(a(t))£- v / 2Re(a(t))^ 


(5.139) 


On the other hand, as a check, we can also arrive at this result by starting from ex¬ 
pansion (5.125), where we interpret the sum over n as the power series for the exponential 
map, i.e. 


1 I ■ 9 iut 

t (t,x ) = e — 21 "! e ~ 2 


E 

n =0 


(a(f)a t ) n _ 1 i a |2 

————= e 21^1 


Til 


jut / , \ + 

e -^- e a(t)at 


(5.140) 
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To show that this result is the same as Eq. (5.139) we simply calculate 


4t)A-aHt)a Mx) Eq. J2.72) ^ a t e -«*( t )a J aa * [<J , a] ^ = 

e -|l«l 2 gaWot _ a *^ + ...) = 

e -|l«| 2 gaWat s (5.141) 


where we once again used the fact, that the action of the annihilation operator on the 
ground state vanishes. Inserting the result into Eq. (5.140) we get 


ILdt . . j . , . 

t (t,x) = e -^ e a{t)a a ■ 


(5.142) 


Explicit calculation of the wave function: Our aim is now to explicitly compute 

the wave function from expression (5.138). We first state the result and derive it in the 
following 


— 


1 


icot 

e~ 2 e 


y/2 Hm(a(t)) £ — \/2Re(a(t)) 


d 

e 2 


7T 


1_ e “ TT 2 a(f) £ - y - Re(a(i)) a(i) 


(5.143) 


7T 


Proof: Let us now perform the explicit calculation. That means, we apply the 

differential operator 7 e v = e ^ llm OO)i V 2 Re(a(t)) ^ ^he ground state of the harmonic 

S? 


oscillator ipo(£) = 


\Jxo pv 

the operator expansion, as sketched in Eq. (5.144) 


e 2 , but we will only consider the linear and quadratic term of 


e v e~ i 2 = (1 + V + i V 2 + 


e 2 , 


(5.144) 


linear term: 


C / \ e 2 

De”2 = ( \[2i Im(o:(f)) £ — \/2 Re(o:(f) e~ 

= ^\/2 i Im(a(t)) £ + \/2 Re(o:(f) e _ 2 

= \/2£ (Re(a(t)) + i!m(o:(f))) e~ 2 


«(*) 


= \/2 ct(£) £ e 2 , 


(5.145) 


' The operator T >, which is a notation to shorten the cumbersome calculations, is not to be confused 
with the displacement operator D. 
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quadratic term: 


V 2 e~ 2 


j- C 

V2a(t)V£e ~ 2 (5.146) 

\Pl a(t) (V2i Im(o:(t)) £ — V2Re(a(t ))£ e~ 2 

(2 a(t) i Im(o:(t)) £ 2 — 2 a(t) Re(a(t)) + 2 a(t) Re(o:(t)) £ 2 ) e 2 

£ 

(2a(t) (Re(a(t)) + ilm(a(t))) £ 2 — 2 a(t) Re(a(t))) e 2 . 

"-——.-' 

a(t) 


Inserting the results of Eq. (5.145) and Eq. (5.146) into the expansion (5.144) we get 
e v e~^ = (1 + V2a(t)£ - a(t)Re(a(t)) + ...) e “V, (5.147) 

v -v'-' 2! 

linear term 

where we can already identify the linear term and the Erst part of the quadratic term in 
the end result of the expansion. Finally, rewriting the whole expansion as exponential 

function and including the factor e 2 , we arrive at the proposed result of Eq. (5.143). 


Probability distribution of the wave packet: As last computation in this section 

we want to study the probability density of the wave function of the coherent state, thus 
we take the modulus squared of expression (5.143) 

x)\ 2 = 4>* a (f> a {t,x) 

2Re(a(t)) 2Re(cc(t)) 

_ 1 e -£ 2 + V2 ( a* + a) £-Re(a(t)) (a* + a ) 

Xo yCl 

— ^ c ~(g~ y 2 Re(a(t ))) 2 _ 1 c ~(g~ v^lal cos (ujt-<p )) 2 (5148) 

x oy /n x 0 

Recalling the result for the expectation value (the mean value) of the position (Eq. (5.134)) 
and our choice of £ = — we finally hnd the probability density 

Xq 

\0a(t,x)\ 2 = —exp - -f- — ^ , (5.149) 

\ Xq j 

which is a Gaussian distribution with constant width. 


Result: The coherent states are oscillating Gaussian wave packets with constant 

width in a harmonic oscillator potential, i.e., the wave packet of the coherent state is 
not spreading (because all terms in the expansion are in phase). It is a wave packet 
with minimal uncertainty. These properties make the coherent states the closest quantum 
mechanical analogue to the free classical single mode field. For an illustration see Fig. 5.3. 
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Figure 5.3: Coherent state: The probability density of the coherent state is a Gaussian 
distribution, whose center oscillates in a harmonic oscillator potential. Due to its being a 
superposition of harmonic oscillator states, the coherent state energy is not restricted to 
the energy levels hu(n + |) but can have any value (greater than the zero point energy). 





